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Abstract. We review our recent achievements in the construction of microscopic mass tables based on the 
Hartree-Fock-Bogoliubov method with Skyrme effective interactions. In the latest of our series of HFB- 
mass models, we have obtained our best fit ever to essentially all the available mass data, by treating 
the pairing more realistically than in any of our earlier models. The rms deviation on the 2149 measured 
masses of nuclei with N and Z > 8 has been reduced for the first time in a mean field approach to 0.581 
MeV. With the additional constraint on the neutron-matter equation of state, this new force is thus very 
well-suited for the study of neutron-rich nuclei and for the description of astrophysical environments like 
supernova cores and neutron-star crusts. 

PACS. 21.10.Dr Binding energies and masses - 21.30.Fe Forces in hadronic systems and effective inter- 
actions - 21.60.Jz Hartree-Fock and random-phase approximations - 26.60.Gj Neutron star crust 



1 Introduction 



latter constraint is not only of prime importance for re- 
liable investigations of a possible superfluid phase in the 
inner crust of neutron stars, but also turns out to signifi- 
cantly improve the accurary of the mass fit. 



2 The Hartree-Fock-Bogoliubov mass 
formulas 



force of the form 



„Sky 



Nuclear astrophysics applications require the knowledge 
\ of various nuclear properties (nuclear masses, nuclear level 
densities (NLD), optical potentials, 7-ray strength func- 
tions, etc.) that cannot be measured experimentally in the 
foreseeable future. In order to make reliable extrapolations 
of these quantities far from the domain covered by exper- 
imental data, we have developed a series of nuclear-mass 
. models based on the Hartree-Fock-Bogoliubov (HFB) method^ Hartree-Fock-Bogoliubov (HFB) mass models HFB- 
; with Skyrme and contact-pairing forces, together with phe- 16 and HFB- 17 are based on the conventional Skyrme 
nomenological Wigner terms and correction terms for the 
spurious collective energy. The model parameters are fit- 
ted to essentially all the available atomic mass data, re- 
quiring that the model reproduce several properties of 
uniform asymmetric nuclear matter determined by mi- 
croscopic calculations with realistic nucleon-nuclcon po- 
tentials. With these nuclear-matter constraints, our mod- 
els can be reliably applied to study astrophysical environ- 
ments such as supernova cores and neutron star crusts [1] . 
In particular, model HFB-9 [2] and all later models con- 
strained the underlying Skyrme force to fit the equation 
of state of neutron matter, as calculated by Friedman 
and Pandharipande [3] for realistic two- and three-nucleon 
forces. 

In this paper, we present our latest models, HFB-16 [4] 
and HFB- 17 L S\ 7 in which we have imposed the additional 
constraint of reproducing as a function of density the 1 S'o 
pairing gaps of uniform asymmetric nuclear matter. This 
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where r tj = r i -r J ,r = (r l +r j )/2, p y = -ih(Vi- Vj)/2 
is the relative momentum, and P a is the two-body spin- 
exchange operator. The contact pairing force acts only 
between nucleons of the same charge state q (q = n or p 
for neutron or proton, respectively) and is given by 
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where v n q [p ni p p ] is a functional of the nucleon densities. 

To the HFB energy calculated for the Skyrme force (p} 
and the pairing force ([2]) are added two phenomenological 
corrections: (i) the Wigner energy [6j[T] 



Ew = Vw exp 
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and (ii) the rotational and vibrational spurious collective 
energy 



E co \\ — 
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in which E^ nk denotes the cranking-model value of the 
rotational correction [5] and /3 2 the quadrupole deforma- 
tion, while all other parameters are free fitting parameters. 
The correction term, Eq , differs from that used in our 
previous mass models. While the differences are small for 
large deformations, the collective energy given by Eq |J4]) 
now vanishes for spherical nuclei. 



3 Microscopically derived effective pairing 
force 

In all of our mass models HFB-1 to HFB- 15, the density 
dependence of the effective pairing strength v 77 q [p n , p p ] 
acting between nucleons of the same charge q (q — n or 
p for neutrons or protons, respectively) was parametrized 
by an expression of the usual following form 



■>* q \p n ,p P ]=V 1[ Jl-ri. (£) } 



(5) 



where p = p n + p p is the total density, po is the satura- 
tion density of symmetric nuclear matter, whereas V nq is 
a free parameter that were determined by the global mass 
fit. The surface parameters , r\ s and a, were either set to 
zero assuming a purely volume pairing force or taken from 
the work of Ref. [9j . Even when constraining the pairing 
strength V 7rq to a relatively low value in order to con- 
form the nuclear pairing with odd-even mass differences, 
we found that our earlier mass models predict unrealis- 
tic pairing gaps in uniform nuclear-matter as compared to 
microscopic calculations while yielding comparably good 
mass fits (hence comparably good pairing gaps in finite 
nuclei). 

In order to improve the treatment of pairing, we have 
recently developed a new series of HFB mass models by 
imposing the additional constraint of reproducing as a 
function of density the 1 Sq pairing gaps of uniform asym- 
metric nuclear matter. Fitting all the model parameters 
while imposing this constraint would be an extremely oner- 
ous numerical task when using the ansatz (|5|). Instead of 



postulating a density dependence as in Eq. ([5]), we have 
thus determined the strength of the effective pairing force 
at each neutron and proton density by solving the HFB 
equations in uniform matter and requiring that the re- 
sulting gap reproduce exactly the microscopic pairing gap 
calculated with realistic forces at that neutron and proton 
density. In this way, the pairing strength for the nucleon 
species q is given by 



3/2 
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where A q {p ni p p ) is the corresponding pairing gap in uni- 
form matter, M*(p n , p p ) is the effective nucleon mass and 
£a is the pairing cutoff. The chemical potential p q is ap- 
proximated by 



2 1.2 



Pq = 



Fq 
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where k-Fq = (Sir 2 Pq) 1 ^ is the nucleon Fermi wave num- 
ber. 

Given the discrepancies between the various micro- 
scopic many-body calculations of the 1 Sa pairing gaps in 
uniform nuclear matter, we have considered two different 
cases. In model HFB-16, we have taken the microscopic 
pairing gap calculated at the lowest BCS level using the 
Argonne vl4 potential 10J and following Duguet [llj . we 
have assumed that the pairing force for nucleons q de- 
pends only on the density p q . We have dropped this as- 
sumption in our latest model HFB- 17 and we have consid- 
ered the recent Brueckner calculations of Ref. [12] which 
include the effect on the interaction of medium polariza- 
tion. These calculations were performed using the Argonne 
vl8 potential both with and without self-energy correc- 
tions. Adopting a microscopic pairing gap calculated with 
self-energy corrections is quite challenging since ideally for 
consistency the Skyrme effective mass M* should be fit- 
ted to the corresponding microscopic effective mass. How- 
ever this seems to be impossible within the framework of 
the conventional Skyrme forces used here. For this reason, 
we have used the microscopic pairing gap of Ref |12j in- 
cluding only medium-polarization effects and we have set 
M* = M in Eqs. © and 0. 

Ref. (12] calculates pairing gaps only for symmetric nu- 
clear matter, A$m(p = Pn+Pp), and pure neutron matter, 
Anm{Pu)- Since we need the pairing gaps for arbitrary 
asymmetry we adopted the interpolation ansatz 

A q (p n , Pp ) = A SM (p)(l - M) ± A NM ( Pq ) Vj, (8) 

where i] = (p n — p p ) / p and the upper (lower) sign is to 
be taken for q = n(p); we have also assumed charge sym- 
metry, i.e., A n (p n ,p p ) = A p (p p ,p n ). This expression en- 
sures that for symmetric nuclear matter, A q (p/2, p/2) = 
Asm(p) and for neutron matter A n (p, 0) = Anm(p) and 
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Table 1. Parameters of the analytical fit for the microscopic 
reference pairing gaps (the unit of length is fermi and the unit 
of energy is MeV) . 
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k 2 


k 3 


km 


SM 


f33.779 


0.943146 
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2.11577 


1.51 


NM 


14.9003 


1.18847 


1.51854 


0.639489 


1.52 
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Fig. 1. 1 So neutron pairing gap A in infinite neutron matter as 
a function of the density p. The dashed line represents the gap 
obtained in the BCS approximation from Ref. |10) and used 
in model HFB-16, while the solid line is the gap calculated 
in Ref. [12] including the effect on the interaction of medium 
polarization and used in model HFB-17. 



A p (p, 0) = 0. For convenience, we have used the essentially 
exact analytical representations 



Asm(p) = 0(k m - hp) A 
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A NM (p n ) = 0(k m - k Fn ) A 



k 2 



°~ky n + k 2 (k Fn - k 2 ) 2 + k 2 ' 

(10) 

where k F = (37r 2 p/2) 1 / 3 , 9 is the Heaviside unit-step func- 
tion, and the associated parameters are given in Table [TJ 
The different pairing gaps are shown in Figs. [1] and [D 
For odd-nucleon nuclei we make use of the "equal-filling" 
approximation |13j . In order to take phenomenologically 
into account the contribution to pairing of the time-odd 
fields [14] , we multiply the pairing strength by a parity fac- 
tor (+ for an even number of nucleons, — otherwise) 
that we allow to be different for neutrons and protons due 
to Coulomb forces and possible charge symmetry breaking 
effects. By definition we set /+ = 1. 

The corresponding effective pairing force, as given by 
Eq. @, is shown in Figs. [3] and H for the HFB-17 mass 
model. 




"0 0.04 0.08 0.12 0.16 0.2 0.24 0.28 
P [fm -3 ] 

Fig. 2. 1 5'o neutron pairing gap A in infinite symmetric nu- 
clear matter as a function of the density p. The dashed line 
represents the gap obtained in the BCS approximation from 
Ref- ED] an d used in model HFB-16, while the solid line is the 
gap calculated in Ref. [12] including the effect on the interac- 
tion of medium polarization and used in model HFB-17. 
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Fig. 3. Isospin dependence of the effective neutron pairing 
strength v nn of model HFB-17 for different densities p (indi- 
cated above each curve in fm -3 ). 
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Fig. 4. Effective neutron pairing strength v^ n of model HFB- 
17 as a function of density p for different isospin asymmetries 
r\ by step of 0.2. 



4 Results of the mass fit 



As explained in the previous sections, the HFB-17 mass 
model [S] represents our latest attempt to provide a uni- 
versal force capable of predicting globally static nuclear 
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Fig. 5. Differences between experimental and calculated 
masses as a function of the neutron number N for the HFB-17 
mass model. 



structure properties. The 10 parameters of the Skyrme 
force along with the 4 parameters of the pairing force (in- 
cluding the pairing cutoff) and the 9 parameters of the 
Wigner and collective corrections have been fitted on the 
2149 measured masses of nuclei with N and Z > 8 given 
in the 2003 AME [15] . The deviations between the experi- 
mental data and the HFB-17 predictions are shown graph- 
ically in fig. El The rms and mean (data - theory) values of 
these deviations are 0.581 MeV and -0.019 MeV, respec- 
tively. HFB-17 is not only the most accurate mass model 
ever achieved within the mean-field framework, it also sat- 
isfies extra physical constraints that make it more suitable 
for astrophysics applications. In particular, the calculated 
quadrupole moments, charge radii, charge-density distri- 
butions and spins are also found to be in good agreement 
with experiment. As done in our previous mass fits, when 
determining the Skyrme force parameters, a special atten- 
tion has been paid to properly describe the properties (not 
only the pairing) of infinite nuclear and neutron matter de- 
termined from realistic calculations [3,16,17]. The excel- 
lent fit of HFB-17 to the neutron matter curve of Ref. [3] 
was achieved simply by imposing a nuclear-matter sym- 
metry coefficient of J = 30 MeV. The potential energy 
per particle determined with the BSkl7 Skyrme force was 
also found to be in fair agreement with the recent calcula- 
tions of Ref. [16] in each of the four two-body spin-isospin 
(S, T) channels and for all densities. Finally, in contrast 
to most of the traditional Skyrme forces, our BSkl7 (and 
previous forces) is consistent with the isovector splitting of 
the effective mass deduced from measurements of isovec- 
tor giant resonances and confirmed in several many-body 
calculations with realistic forces [17] . More details on these 
comparisons can be found in Ref. 

In Fig. [5] we compare the HFB-17 predictions with 
those of the finite range droplet model (FRDM) [18] for 
all 8 < Z < 110 nuclei lying between the proton and neu- 
tron drip lines. Significant differences (up to 25 MeV) are 
found, especially for heavy neutron-rich nuclei. Such dif- 



30 
25 
20 
> 15 

2 10 

< 


-5 
-10 



M(FRDM)-M(HFB-17) 





o 



40 



80 



120 

N 



160 



200 



240 



Fig. 6. Differences between the FRDM and HFB-17 mass pre- 
dictions as a function of (left) N for all 8 < Z < 110 nuclei 
lying between the proton and neutron drip lines. 
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Fig. 7. Differences between the HFB-16 and HFB-17 mass 
predictions as a function N for all 8 < Z < 110 nuclei lying 
between the proton and neutron drip lines. 



ferences are much larger than what is obtained between 
our various HFB mass models, as illustrated in Fig. \7\ 

In Figs. [5] and [HI we show the neutron-shell gaps, de- 
fined by 



A n (N 0) Z) = S 2n (N , Z) - S 2n (N + 2, Z) 



(11) 



as a function of Z for the magic numbers Nq — 50, 82, 126 
and 184 for the HFB-17 model (SVn is the two-neutron sep- 
aration energy). Likewise in Fig. [10] we show the proton- 
shell gaps as a function of N for Zq = 50 and 82. Here 
we consider just the comparison with the data. Given the 
overall close agreement between experimental and theo- 
retical shell gaps, the differences found for the Nq = 126 
neutron-shell gaps clearly emerge as unsatisfactory, espe- 
cially in the vicinity of doubly-magic 208 Pb nuclei. It is 
even more disturbing that the Zq = 82 proton-shell gaps 
agree rather well. The double magicity around 208 Pb re- 
mains unresolved in all mean field calculations, regardless 
of the type of force or its parametrization. 
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Fig. 8. No — 50 and iVo = 82 shell gaps as function of Z for 
the mass models HFB-17. 
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Fig. 9. No = 126 and No = 184 shell gaps as function of Z for 
the mass models HFB-17. 
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Fig. 10. Zo — 50 and Zo = 82 shell gaps as function of N for 
the mass models HFB-17. 
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Fig. 11. Proton number Z and neutron number N of nuclei in 
the outer crust of neutron stars for the HFB-16 and HFB-17 
mass models. 



5 Astrophysical applications 

Following the classical paper of Baym, Pethick and Suther- 
land [19], we have calculated the structure of the outer 
crust of neutron stars for the HFB-16 and HFB-17 mass 
models. The only microscopic inputs are the values of 
the atomic masses. We have used experimental data when 
available. The results are shown in Fig. QT] The compo- 
sition of the nuclei in the outer crust is essentially the 
same for the two mass models HFB-16 and HFB-17. Both 
models predict that at densities above /?nd — 3.7 x 10 11 
gem -3 , neutrons start to drip out of nuclei. The corre- 
sponding layer of the outer crust is composed of nuclei 
with Z = 38 and N = 82. 

The structure of the inner crust of neutron stars, at 
densities above pnd, is not simply determined by atomic 
masses due to the presence of the neutron gas. The free 
neutrons modify the properties of the "nuclei" by i) ex- 
erting a pressure on them and ii) by reducing their sur- 
face tension 1 . In order to make reliable predictions of 
the composition of the neutron star crust beyond neu- 
tron drip, both the nucleons bound inside clusters and the 
free neutrons have to be described consistently. For this 
purpose, we have applied the Extended Thomas-Fermi 
method up to the 4th order (see Ref. [20] for details) 
with the Skyrme effective forces BSkl6 and BSkl7 under- 
lying the HFB-16 and HFB-17 mass models respectively. 
As can be seen in Fig. [12l the two models yield very simi- 
lar results. The properties of the neutron star crust (elas- 
tic constants, electric and thermal conductivities, etc.) are 
strongly dependent on the proton number Z of the nuclear 
clusters. As shown in Fig. 1131 the differences in Z between 
the two crust models do not exceed 2 units. 



6 Conclusions 

Despite all the restrictions imposed on the model parame- 
ter space, we have obtained the most accurate mass table 
ever achieved within the mean field framework, the rms 
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Fig. 12. Proton number Z and mass number A of "clusters" 
in the inner crust of neutron stars for the HFB-16 and HFB- 
17 mass models. Note that A includes not only the number of 
bound nucleons but also the mean number of free neutrons per 
lattice site. 
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Fig. 13. Proton number Z of nuclei in the inner crust of neu- 
tron stars for the HFB-16 and HFB-17 mass models. 
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deviation falling to 0.581 MeV. Given also the constraint 
imposed on the Skyrme force by microscopic calculations 
of nuclear and neutron matter, this new model is partic- 
ularly well adapted to astrophysical applications involv- 
ing a neutron-rich environment, such as the elucidation 
of the r-process of nucleosynthesis, and the description 
of supernova cores and neutron-star crusts. The HFB-17 
mass table is made available to the scientific community 
at |http://www-astro.ulb.ac.be| 



References 

1. N. Chamel and P. Haensel, "Physics of Neu- 
tron Star Crusts", Living Rev. Relativity 11, 
(2008), 10. URL (cited on 15 D ecember 2008): 
http: / /www. livingreviews.org/lrr-2008-10 

2. S. Goriely, M. Samyn, J. M. Pearson, and M. Onsi, Nucl. 
Phys. A750 (2005) 425. 

3. B. Friedman and V. R. Pandharipande, Nucl. Phys. A361, 
502 (1981). 

4. N. Chamel, S. Goriely, and J. M. Pearson, Nucl. Phys. 
A812, 72 (2008). 



